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ABSTRACT 
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CO Context. Recently, much attention has been paid to the magnetic couphng (MC) process, which is supported 

(N ■ 

■ by very high emissivity indexes observed in Seyfert 1 galaxy MCG-6-30-15 and GBHC XTE J1650-500. 
' But the rotational energy transferred from a black hole is simply assumed to be radiated away from the 
' surrounding accretion disk in black-body spectrum, which is obviously not consistent with the observed 
Oh hard power-law X-ray spectra. 

O Aims. We intend to introduce corona into the MC model to make it more compatible with the observations. 
' Methods. We describe the model and the procedure of a simplified Monte Carlo simulation, compare the 
output spectra in the cases with and without the MC effects, and discuss the influences of three parameters 



involved in the MC process on the output spectra. 

Results. It is shown that the MC process augments radiation fluxes in the UV or X-ray band. The emergent 
, spectrum is affected by the BH spin and magnetic field strength at the BH horizon, while it is almost 
unaffected by the radial profile of the magnetic field at the disk. 

Conclusions. Introducing corona into the MC model will improve the fitting of the output spectra from 
AGNs and GBHCs. 
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1. INTRODUCTION 

Tenuous hot plasma (corona) is commonly believed to exist in the inner region of the ac- 
cretion flow, which is produced and heated by magnetic activities, such as reconnection and 
flare etc. (e.g. {Liang & Price 1977) IGaleev et al. 1979V The influences of the disk corona on 
the emergent spectrum from active galactic nuclei (AGNs) and galactic black-hole candidates 
(GBHCs) have been discussed by many authors (e.g. Haardt & Maraschi 1991, hereafter lHM9 II 
IHaardt & Maraschi 19931 {Field & Rogers 1993| lEsin 19991 {Zhang et al. 2000| ILiu et al. 20021 . 
Corona explains the power-law X-ray spectra very well, and reprocessing of the coronal X-rays 
Send offprint requests to: D. X. Wang, e-mail: dxwangOhust . edu . cn 
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by cold disk gives rise to the observed emission lines naturally, in which iron Ka fluorescence 
line provides us a diagnostic of the geometry of the accretion flow and the property of the space- 
time around the black hole (BH). However, only gravitational energy of the accreting matter is 
invoked to heat the disk and its corona in previous works. 

Recently, much attention has been paid to the magnetic coupling (MC) process in which en- 
ergy is transferred from a rotating BH to its surrounding disk JBlandford 1999llvan Putten 19991 
|Li & Paczynski 2000| Li l2002al Wang et al. l200l Wang et al. 2003a, b, hereafter IW03al and 
IWU3bl respectively), which is in fact a variant of the Blandford-Znajek (BZ) process. The exis- 
tence of the MC process has been supported by very high emissivity indexes observed in Seyfert 
1 galaxy MCG-6-30-15 and the GBHC XTE J1650-500, which cannot be explained by flie stan- 
dard accretion disk (Wilms et al. 120011 iMiTler et al. 2002llLi 2002bt . 

However, in the previous MC models all the energies that deposit on the disk, both the grav- 
itational energy of the accreting matter and the rotational energy transferred from the BH, are 
simply assumed to be radiated away in black-body spectrum, which is obviously not consistent 
with the observed hard power-law X-ray spectra. In this paper we intend to introduce corona 
into the MC model and discuss the spectrum affected by the MC process. It is shown that the 
MC process augments the radiation fluxes, while corona makes MC model more reasonable and 
more consistent with the observations than the previous MC model. 

The paper is arranged as follows: Sect. 2 gives the description of our model. Sect. 3 shows 
the procedure and results of our simplified Monte Carlo simulation and in Sect.4 we discuss the 
applications of our model to observations. Throughout this paper the geometric units G = c = 1 
are used. 



2. DESCRIPTION OF OUR MODEL 

The poloidal profile of the large-scale magnetic field and the geometry of the corona are illus- 
trated in Fig.[n The central BH is fast rotating and surrounded by axisymmetric magnetosphere. 
It is assumed that the disk is geometrically thin and optically thick, while the corona is geo- 
metrically thick and optically thin. The poloidal magnetic field is assumed to be constant on the 
horizon and it varies as a power-law with the radial coordinate on the disk, i.e., 

Bh = const, Bd oc r'". (1) 

In Fig. 1 the angle 9m indicates the angular boundary between the open and closed field lines 
at the BH horizon, which can be determined by the mapping relation derived in lW03bl The angle 
Ol is the lower angular boundary of the closed field lines, and it is supposed to be less than O.Stt 
to avoid the singularity of the closed field lines at the equatorial plane. The magnetic field powers 
the jet via the BZ process and transfers energy between the BH and the disk via the MC process. 

Both large- and small-scale magnetic fields are involved in our model though the latter is not 
shown in Fig.^ The large-scale magnetic field plays a key role in transferring energy and angular 
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momentum from a fast-rotating BH to its surrounding disk or remote astrophysical load in the 
MC and BZ processes, respectively. While the small-scale magnetic field energizes the corona 
above the disk by virtue of some physical processes, such as magnetic reconnection and buoyancy 
JHM91> . Moreover, the large-scale field can be produced from the small-scale field created by 
dynamo processes (e.g. [Tout & Pringle 1996| l. Livio et al. J1999t argued that the strength of 
the large-scale field threading a black hole is very weak if the field is created in the thin disk, 
because the ratio of the large-scale field to the small-scale field is ~ HjR <s; 1, where H is the 
disk thickness and R is the radial size of the field. This greatly confines the BZ and MC powers. 
Since the thickness of the corona is about the radius of the disk, the large-scale magnetic fields 
created by dynamo processes in the corona are significantly stronger than those in the thin disk 
JCao 2004> . and the strong BZ and MC powers are maintained. For simplicity, we suppose that 
these two kinds of magnetic fields exist independently. 

The energy deposited in the corona and disk consists of two parts: (i) the gravitational energy 
of the accreting matter, and (ii) the rotational energy of the BH transferred in the MC process. 
Some authors (e.g. lHM91IIMerloni & Fabian 2002t have mentioned physical processes in which 
the gravitational energy converts to the thermal energies of the corona and disk. However, nobody 
has proposed any mechanism to describe how the rotational energy of the BH transferred in the 
MC process dissipates to date. In this paper, we deal with this issue based on the following two 
considerations. 

(i) Based on BH magnetosphere theory JMacdonald & Thorne 19821 IThorne et al. 19861 Li 
l2002allW03a> . the disk can be regarded as an electromotive force (EMF) in an equivalent circuit 
due to its rotation, which is opposite to the BH EMF. If the angular velocity of the BH is greater 
than that of the disk, Poynting flux would flow from the BH to the disk JW03bt . Although the 
disk is commonly supposed to be a perfect conductor, its interior resistance may not be zero. Part 
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of the MC power may heats the disk, just as Poynting flux gets into the resistance and heats it in 
a ckcuit; 

(ii) The flux of the angular momentum transferred magnetically from the BH to the disk 
decreases rapidly with disk radius (see Eq.(|3} below), and it results in the differential rotation of 
the accreted matter. Thus the electromagnetic energy converts to thermal dissipation, just as the 
gravitational energy of accreting matter converts to thermal dissipation in disk accretion. 

Detailed discussions of these physical processes are beyond the scope of this paper, and we 
assume that both the gravitational energy of the disk and the rotational energy of the BH can 
dissipate into thermal energy of the corona and disk. 

The scenario of the emergent spectrum is described as follows. Thermal soft seed photons 
are emitted from the disk, and when they cross the rarefied hot corona, some may go through the 
corona directly without being scattered, some may escape after one or several times of scatter- 
ing, and others may be down-scattered to the disk by the corona. The fate of the X-ray photon 
that returns to the disk depends critically upon its energy E: (a) if £ > lOOkeV, the photon 
loses its energy by Compton recoil, (b) if E < a few keV, the photon is absorbed by ionization 
or free-free absorption, and (c) it may be both Compton scattered and bound-free absorbed for 
a /ewkeV < E < lOOkeV ( [Field & Rogers 1993^ . Bound-free absorption or ionization results 
in a vacancy in the inner shell of an atom, which leads to a fluorescence photon or the ejection 
of an Auger electron (for details see George & Fabian 1991, hereafter IGF91> . Unscattered pho- 
tons form the component of the UV/soft X-ray band, photons escaping the corona after inverse 
Compton scatterings constitute the power-law X-ray spectrum, while reprocessed photons by the 
disk produce fluorescence lines and a reflection hump. 

Here we neglect cyclo-synchrotron radiation since it is less important than Compton ra- 
diation in the slab corona ( IDi Matteo et al. 19971 IGhisellini et al. 1998I I. We also neglect the 
bremsstrahlung radiation, assuming that the corona is tenuous enough. 

From the above scenario it can be found that part of the energy dissipated in the corona 
is conveyed to the disk via hard X-rays, i.e., the disk is heated by three kinds of energies: (1) 
gravitational energy of the accreting matter, (2) rotational energy of the BH transferred to the disk 
and (3) the absorbed energy of X-rays coming from the corona. Since the energy of the X-rays 
comes eventually from the gravitational and rotational energies, we suppose the fraction of the 
total energy dissipate in the disk to be rj, which includes the contribution of the downward X-rays. 
As shown in lHM9ll the soft emission is derived almost only from absorption and reprocessing 
of the high-energy flux impinging on the disk. Roughly half of the Comptonized photons in 
the corona irradiate and heat the disk, and the rest half radiate away. So we take rj - 0.5 in 
calculations. 

Additionally, for simplicity of simulation, the disk is assumed to be "cold", with hydrogen 
and helium being fully ionized, but all other elements being neutral or weakly ionized. Moreover, 
the energetic and spatial distributions of free electrons in the corona is postulated to be thermal 
and homogeneous. 
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3. SIMULATION OF DISK SPECTRUM 

3.1. Multicolor Black-Body Spectrum of the Disk Affected by the MC Process 

Based on the conservation laws of mass, angular momentum and energy, the local radiation flux 
from the disk around a rotating BH is given by Page & Thorne <1974> and Li ( I2002a> . 

F^^f + (-^)(E^ - ^loL^r' r (E^ - nr>L^)Hrdr, (2) 
4Kr rar J^_^^^ 

where Md, ^d, fms, and are respectively the accretion rate, the angular velocity of the 
disk, the radius of the marginally stable orbit, the specific energy and angular momentum at the 
disk radius r. The function / is related to the radiation flux due to disk accretion, and H is the 
flux of angular momentum transferred between the BH and the disk (IW03a> . 

, 1 < 4^ < ^out 

H(a,-^,n)/Ho^\ , (3) 

^ > ^out 




In [2 csc2 0-(l-q)]^| 1 - 2x-^r' + alx-n,ir' ' 
where Ho = (B^)^M, /3 is the ratio of the angular velocity of the disk to that of the BH, i.e. /3 = 



Q.d/^h, ^ = r/rms is a radial parameter of the disk defined in terms of r,„s , Sind Xms = y/fms/M. 
The relation between the magnetic field and the accretion rate is given by Moderski et al. 



J1997> based on the balance between the pressure of the magnetic field on the horizon and the 
ram pressure of the innermost parts of an accretion flow, i.e., 

(BP^f/m = Pram ~ pc^ ~ Mo/(47r4), (5) 

where r// is the radius of the horizon. From Eq.(|5jl we can define = (B^)^ = 
2Mzj/[M2(i + y/l-alf]. 

For the presence of corona the dissipated power in the unit area of the disk, Fj, is related to 
F(r) by 

Fd(r) = r]F(r), (6) 

where F{r) is the local radiation flux expressed by Eq.(|2|i. According to Stefan-Boltzmann law 
we have the local effective temperature on the disk expressed by 

Tdr) = (Fd/(TsBy'\ (7) 

where ctsb is the Stefan-Boltzmann constant. The local radiation spectrum is defined by Planck 
function: 
2h 

c2 exp [hv/kBTd(r)] - 1 
Thus the multicolor black-body spectrum of the disk is 

L, = B„(r)2nrdr. (9) 
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3.2. Monte Carlo Simulation 

If the disk is covered by corona, we need to resolve the radiative transfer in the hot 
corona to get the emergent spectrum. Radiative transfer has been computed in different ways. 
One approach is to solve the radiative transfer equation either numerically or analytically 
(e.g. [Sunyaev & Titarchuk 1980| IPoutanen & Svensson 1996> . Another approach is the Monte 
Carlo simulation (e.g. Pozdnyakov et al. 1983, hereafter IPS 3 1 IGorecki & Wilczewski 19841 
IGF9Ul5temetal. 19951 Hua 1997, hereafter IW71 Yao et The emergent spectrum 

can be also obtained by some approximate semi-analytical formulae (e.g. IZdziarski 19861 
IHua & Titarchuk 1995t . In this paper we calculate the spectrum by using the Monte Carlo simu- 
lation. 

Compared with H97, photon absorption in the disk is taken into account in our model, which 
makes the simulation slightly more complex. The steps of the simulation are: (i) sample a seed 
photon including its position, energy and direction; (ii) draw a value for its free path and test 
whether it can leave the disk or corona; (iii) simulate the interaction of the photon with the 
medium; (iv) repeat steps (ii), (iii) till the photon leaves the system of the corona and disk. 

3.2.1 . Sampling the seed photons 

The probability density of the seed photon can be written as 

p(r, v) = Inrdr ■ F,(r)/L, (10) 

where Fy,(r) is the flux density at r, L = J^'°" L(r)dr is the total luminosity of the disk, and 
L{r)dr - 27irFd{r)dr is the luminosity of the ring r ~ r + dr. Obviously, p(r, v) can be expressed 
as follows, 

L(r)dr 2nr-F,.{r)dr L{r)dr F,.{r) ,,,, 

p{r, v) = — — = — — ■ — — = p(r)B,,(r) (11) 

L L(r)dr L Fdin 

where p{r) is the probability of a photon emitted in the ring r ~ r + dr, which can be sampled by 
tabulation. The Planck function B^,(r) can be sampled in the way described in lP83l 



3.2.2. Drawing the free path 

For the photon scattered in the corona, the probability it travels at least an optical depth r is 
e so the optical depth that the photon travels between the i-th and (i-i-l)-th scatterings can 
be drawn with r, = -In A, where < A < 1 is a random number corresponding to a random 
event. Then the free path of the photon can be drawn with ^ — ■ where rig, cr and 
CTj - 6.65 X 10"^'*cm^ are the number density of the electrons, cross-sections of scattering and 
Thomson Scattering, respectively. The parameters Hg and are the vertical height and optical 
depth of the corona, respectively. Because the corona is very hot, electrons in it are relativistic, 
which means that cr depends not only on the energy of the photon but also on the energy and 



R.-Y Ma et al.: Title Suppressed Due to Excessive Length 



7 



direction of the electron. So the cross-section averaged over the distribution of the electrons in 
IH97l is used to draw the free path of the photon. 

Since the disk is assumed to be "cold", with hydrogen and hehum being fuUy ionized, and 
all other elements being neutral or weakly ionized, hard X-ray photons irradiating the disk from 
the corona can be absorbed by the atoms, as well as being scattered by the free electrons. In this 
case, there are two viable ways to draw the free path of the photon: (i) choose the lesser of the 
two free paths that are drawn with (Ta and cr,, where cr„ and cr, are respectively the cross sections 
of absorption and scattering JGF9 1> : (ii) draw the free path with ^/crJcra'+7^s), and determine 
the interaction by drawing a random number A and comparing it with the probability that a free 
path ends with absorption, i.e., ( = ^r+a~' - f ' '^^e photon will be absorbed and otherwise be 
scattered ( |Rybicki & Lightman 1979) . We follow the first way in our simulations. Since electrons 
in the disk is nonrelativistic, ctj can be given by the Klein-Nishina formula: 



where x - Ihv/nieC^ is the energy of the incident photon in unit of electron-rest-energy. The 
cross section of absorption era is taken from Morrison & McCammon (il983l) . 

As the free path is known, the position that the photon arrives before the next interaction 
can be calculated. If the photon is outside the corona-disk system, it will escape away from the 
system and its energy and direction are recorded. But if the photon transfers from the disk to the 
corona or inverse, the point where the trajectory of the photon cross the interface between the 
disk and corona should be regarded as the next initial position of the photon in calculations (see 
lH97l for details). 

3.2.3. Simulating tine interaction 

If the photon is scattered in the corona or disk, we can first sample an electron from thermal 
distribution, and then calculate the energy and direction of the scattered photon, following the 
procedure described in lP83l 

The bound-free absorption of hard X-rays by the atoms in the disk will lead to ionization 
and vacancy, and induce emission of fluorescence lines with the probability called fluorescence 
yield, Y, rather than ejection of Auger electrons. In the simulations, if the photon is absorbed by 
a certain atom or ion, we can draw a random number A (0 < A < 1) and compare it with the 
corresponding fluorescence yield. If A < Y, an emission line is brought out, whose direction can 
be sampled from the isotropic distribution. If A > Y, an Auger electron is sent out, the photon 
vanishes and its trajectory ends. 

Because of the large abundance and cross-section of absorption of iron, only the 6.4keV Fe 
Ka fluorescence line is considered in our simulations. If the energy of the absorbed photon is 
less than the iron K-shell absorption edge, i.e. E < 7.1keV, no emission line is produced, while 
Fe Ka lines emanate with Y - 0.34 for E > 7.1keV. 




(12) 
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3.3. Parameters of Corona and Results of Simulations 

Optically thin thermal Compton scattering spectra are well represented by power-laws 
with exponential cut-ofFs, and the spectral index can be approximately given by (e.g. 
IRybicki & Lightman 197^ 

'--^^ (13, 

InAi 

In Eq.jl3> Ai - 1 + 40 + 160^ is the average photon energy amplification per scattering, where 
= kTclm^c^ with Tc being the coronal temperature. The average scattering probability P for a 
uniform slab is given by |Zdziarski et al. (1994)| as 

exp(-T) / I \ 1 T 
P=l.^(--l)-^.2^:(r), (14) 

where E\ is the exponential integral (IPress et al. 1992> . 

Observations of Seyfert galaxies show that their spectral indexes and the coronal Thomson 
optical depths are about 0.9 and 1, respectively ( IHM91llZdziarski 1998> . and we take a - 0.9, 
T = 1 in simulations. From Ea.J13>. we estimate the temperature of the corona to be Tc ~ 60keV. 

In the following calculations, the radii of the inner and outer edges of the corona and disk are 
taken as r„„ and 20r,„5, respectively. The half heights of the disk and the corona are assumed to 
be Hci-Q.Q2r,ns and Hc-r„,s, respectively. The value of Hd <K r,,,, does not affect the final results. 
The radial range of the disk is taken to be rather small, because the MC power concentrates in 
the innermost region of the disk, and it can be easily verified that at least 80% of the radiation is 
emanated within 20rms- 

Since the typical magnetic field of the supermassive and stellar BHs are respectively about 
10^ G and 10'^ G, we take the middle value, 10^ G, in simulations. A large number (10^) of 
photons are traced to get one spectrum in simulations, and the result of 10^ events costs about 
100 minutes in running the program on a PC with 2.8 GHz Pentium 4 CPU. All the spectra in 
this paper are given in forms of curves rather than histogram to make the curves look smooth. 

The simulation results are given in Fig. 2. The spectrum of the disk-corona system comprises 
three components, the black-body spectrum formed by unscattered photons, the power-law spec- 
trum formed by photons that escape from the corona after several times of inverse Compton scat- 
terings, and the reflected spectrum characterized by iron fluorescence line and reflection hump. 

In order to study the eff'ects of the MC process and corona, the spectra of the disk-corona 
system with and without the MC process and those with and without corona are simulated as 
shown in Fig. 3. 

From Fig. 3 we find: (i) the disk can at most emit soft X-rays with E < a/ewkeV without 
corona, while the spectrum becomes much more complex and much harder when corona exists; 
(ii) the MC process can increase the flux and extend the spectrum to higher energy when the 
corona does not exist, but it can only augment the flux and has no influence on the cutofi" energy 
and the index of the spectrum. 
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Fig. 2. The spectram of the disk-corona system (sohd line) and its components. The black-body, 
power-law and reflected components are shown in dot-dashed, dotted and dashed Unes, respec- 
tively. The parameters a, = 0.998, n = 3 and B = lO^G are taken. 



10000 




hv (keV) 



Fig. 3. Multicolor spectra of the disk and the simulated spectra of the disk-corona system. Thick 
and thin Unes correspond to the spectra with and without the MC process, respectively. Solid 
and dashed lines correspond to the spectra with and without corona, respectively.The parameters 
a. = 0.998, n = 3 and B = lO^G are taken. 

The fact that the cutofi" energy and index of the spectrum are indiff'erent with the MC process 
when corona exists is because the spectral cutofl" hvcm and index a are determined by and r, 
which are independent of the MC process in our model. 

Spectra with difi'erent a„ n and B are given in Fig. 4, from which we find (i) the fluxes 
increase with a, and B no matter the MC process exists or not; (ii) the parameter n has little 
influence on the flux. 



10 
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Fig. 4. Spectra corresponding to different values of the parameters a,, n and B in (a), (b) and 
(c), respectively. Thick and thin Unes correspond to the states with and without the MC process, 
respectively. 



It can be easily explained why the fluxes increase with a, and B no matter the MC process 
exists or not. For pure accretion without the MC process, as a, increases the inner edge of disk 
approaches the BH horizon, and more gravitational energy can be released if the accretion rate is 



R.-Y Ma et al.: Title Suppressed Due to Excessive Length 1 1 

the same. So the fluxes of the disk-corona system increase with a,. According to Eq.Q, stronger 
magnetic field corresponds to a greater accretion rate Mq, and thus higher fluxes. At the presence 
of the MC process, higher fluxes are produced for greater fl» and B since the MC power increases 
with both a, and B. The flux is almost unaffected by n because the MC power is not significantly 
afl'ected by n. 

4. APPLICATION TO ASTROPHYSICS 

In this paper corona is introduced into our model, and it is helpful to interpret the following 
observations. 



4.1. Jet 

The same as our previous model, the BZ and MC processes can coexist in the present model, 
provided that the parameters a* and n are greater than some critical values, and their powers are 
derived by using a modified equivalent circuit and expressed as follows ( IW03b> . 

Pbz/Po = 2al ^ — , O<e<0M, (15) 

Jo ^ " „^-,;„2r 



^0 2-(l-g')sin^ 

j6(l -/?)sin^ft 
2 - (1 - q) sin^ 



P.c/Po^2air^-^^^^1^,0.<0<0., (16) 

J 1 /I „N„.„ . 



where = {B'^^f Ml 

An important feature of our model is the strength of the BZ power relative to the MC power 
can be adjusted by virtue of the angle 0m, which is the angular boundary between the open and 
closed field lines at the BH horizon as shown in Fig.l. As argued in IW03bl 6m is eventually 
determined by the parameters a, and n contained in the mapping relation between the angular 
coordinate on the BH horizon and the radial coordinate on the disk. 

Combining accretion with the MC process, the luminosity of the disk can be expressed as 

L^(1-E,„,)Md + Pmc- (17) 

The first term in Eq.jl7li comes from disk accretion, and the second term from the con- 
tribution of the MC process. The quantity E,„s is the specific energy of the accreting matter 
corresponding to the marginally stable orbit and it reads dNovikov & Thorne 19731 1 



Ml- 



(18) 



Our model can be applied to the observations of some sources, fitting the jet power in radio 
band driven by the BZ process and fitting the bolometric luminosity driven by disk accretion and 
the MC process. 

Here we apply our model to fit the observational data of 3C 273. The BH mass contained is 
estimated as 6.59 x 10'' Mg by reverberation method, the bolometric luminosity L is estimated 
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Fig. 5. The numbers of iron fluorescence photons at different radii and diff'erent incline angles. 
Solid lines, dashed lines and dotted lines correspond to yu = cos / = 0.95, 0.75, 0.35, respectively. 
Thick lines and thin lines coiTespond to the cases with and without the MC effects, respectively. 
The parameters, a» = 0.998, n - 3 and B = lO^G are given. 

to be ~ 10'*^erg ■ s"' iPaltani & Turler 2005> . and the mmimum jet power is given as f ™ = 
6.61 X lO'^'^erg ■ s ' (Celotti et al. ll997> . If the magnetic field at the BH horizon is considered to 
be B''fj ^ lO'^Mg '^^G (Beskin [l997t . we can estimate some parameters of our model as a* ^ 0.76, 
n ^ 6.13 and 6m - 48° by combining Eqs. (I15t - J18l l with the relation (|5}- 

The spectral index, a x 0.77, and luminosity of the X-rays in the range of 3-lOkeV, L3_io ~ 
0.8 1 X 10"*''erg- s"' can be derived from Page et al. J2004t . where some X-ray observational results 
of 3C273 are listed. Combining Ea. (ll3t with the assumption that t - 1, we have Tc - 73keV. 
It can then be obtained by Monte Carlo method that L3_i() ~ 0.72 x 10"*''erg ■ s"', which is a 
little less than the observational data. Since 3C273 is a well known radio-loud quasar with a jet 
showing superluminal motion, the X-rays emitted by synchrotron self-Comptonization of high- 
energy electrons in the jet should be considered. Moreover, the degree of ionization of the disk is 
supposed to be low in our model, which means higher absorption in the range ~ lOkeV. If these 
reasons are considered, the luminosity of the X-rays could be well fitted with our model. 

4.2. Steep Emissivity Index 

Wilms et al. J2001> assumed that the reflected flux has an emissivity profile oc in fitting the 2- 
lOkeV observational data of Seyfert I galaxy MCG-6-30-15, and found a very interesting result. 
The emissivity profile is very steep with /3 - 4.3 - 5.0, which is too high to be explained by the 
standard accretion disk model. 
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In lW03al and lW03EI we explained the very steep emissivity profile by the local radiation flux 
emitted in black-body spectrum by virtue of the MC process. However, the spectra of previous 
models cannot extend to ~ 10 keV for AGNs, which is not consistent with the observations. This 
inconsistency can be removed in the MC model with corona. 

Taking Fe Ka fluorescence as an example, we find in Fig. 5 that the photon numbers produced 
in our model reduce more quickly with the disk radius than those produced in the standard disk 
with corona. It turns out that the very steep emissivity produced in the MC process could not be 
changed at the presence of the corona. 

4.3. QPOs 

Wang et al. J2003ct explained kilohertz quasi-periodic oscillations (QPOs) in BH binaries by 
virtue of rotating hotspots arising from non-axisymmetric magnetic fields. When corona is in- 
troduced into our model, although the rotational energy of the BH is shared by the disk and the 
corona, i.e., the temperature of the hotspot on the disk decreases, the total flux from the region of 
higher magnetic field remains unchanged. So, as the high-magnetic-field region rotates, the fit- 
tings to the frequency of the QPOs are still valid. An interesting phenomenon diflferent with pre- 
vious model is that we can not see a hotter spot in the corona because of the cooling of more soft 
photons from the disk. Additionally, the introduction of corona into our model is helpful to inter- 
pret QPOs which are usually observed in steep power-law state (IMcClintock & Remillard 2003t . 

Although the introduction of the corona into our model is helpful to fitting the observations, 
our model needs to be improved in some aspects. For example, the more reasonable geometry and 
parameters of the corona should be adopted (e.g. lHaardt et al. 19941ILiu et al. 20021 . the cooling 
of synchrotron radiation should be considered, and gravitational effects on the trajectories of 
photons need to be taken into account since it cannot be ignored in the inner region of the disk 
around a rapidly spinning BH. We shall deal with these problems in future works. 
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